Non-Fermi liquid behavior in transport across carbon nanotube quantum dots 
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A low energy-theory for non-linear transport in finite-size single-wall carbon nanotubes, based on 
a microscopic model for the interacting pz electrons and successive bosonization, is presented. Due to 
the multiple degeneracy of the energy spectrum diagonal as well as off-diagonal (coherences) elements 
of the reduced density matrix contribute to the nonlinear transport. A four-electron periodicity with 
a characteristic ratio between adjacent peaks, as well as nonlinear transport features, in quantitative 
agreement with recent experiments, are predicted. 
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Since their recent discovery single-wall carbon nan- 
otubes (SWNTs), of. e.g. m, have attracted a lot of 
experimental and theoretical attention. In particular, as 
suggested in the seminal works 0, Q , due to the pecu- 
liar one-dimensional character of their electronic bands, 
metallic SWNTs are expected to exhibit Luttinger liquid 
behavior at low energies, reflected in power-law depen- 
dence of various quantities and spin-charge separation. 
Later experimental observations have provided a confir- 
mation of the theory 0, 0- As typical of interacting 
electron systems in reduced dimension, SWNTs weakly 
coupled to leads exhibit Coulomb blockade at low tem- 
peratures with characteristic even-odd or four-fold 
periodicity jsl lol llo| . In not only the ground state, but 
also several excited states could be seen in stability dia- 
grams of closed SWNT quantum dots. Such two- fold and 
four-fold character can be qualitatively understood from 
symmetry arguments related to the two-fold band degen- 
eracy of SWNTs and the inclusion of the spin degree of 
freedom. So far, a quantitative description has relied on 
density functional theory calculations {lli] or on a mean 
field description of the Coulomb blockade 
ticular, the position of the spectral lines in the stability 
diagram measured in was found to be in quantitative 
agreement with the predictions in [l^ . However, a mean 
field description may be raoi justified for one-dimensional 
systems. For example, to describe the spectral lines of 
the sample with four- fold periodicity (sample C) in 0, 
a quite peculiar choice of the mean field parameters was 
made, and a quantum dot length three times shorter than 
the measured SWNT length was assumed. Moreover, to 
date no quantitative calculation of the nonlinear current 
across a SWNT dot has been provided. 

In this Letter we investigate spectral as well as dynam- 
ical properties of electrons in metallic SWNT quantum 
dots at low energies. We start from a microscopic descrip- 
tion of metallic SWNTs and include Coulomb interaction 
effects, beyond mean-field, by using bosonization tech- 
niques 0, 13] , yielding the spectrum and eigenfunctions of 
the isolated finite length SWNT. Due to the many- fold 
degeneracies of the spectrum, the current-voltage char- 
acteristics is obtained by solving equations of motion for 



the reduced density matrix (RDM) including off-diagonal 
elements. Analytical results for the conductance are pro- 
vided, which account for the different heights of the con- 
ductance peaks in 0. Moreover, we can quantitatively 
reproduce all the spectral lines seen in sample C in [2j 
by solely using the two ground state addition energies 
provided in that work. The derived level spacing is in 
agreement with the measured SWNT length. 

To start with, we consider the total Hamiltonian 



gate I 



(1) 



where Hq is the interacting SWNT Hamiltonian (cf. 
Eq. |SJ| below) and Hg/j^ describe the isolated metallic 
source and drain contacts as a thermal reservoir of non- 
interacting quasi-particles. Upon absorbing terms pro- 
portional to external source and drain voltages Vs/d 1 they 
read {I = s,d) Hi = Sag ^g.'C^a/CgV/, where ct^; creates 
a quasi-particle with spin a and energy Sq'^i ^ e^— eVg/d 
in lead s/d. The transfer of electrons between the leads 
and the SWNT is taken into account by 



E E /d'^ (r,(f)*t (r)<i>.,(r1 + h-c.) , (2) 

J ^ 



l=s,d <J 



where ^E*^ and ^^^/(r) = <?^^(^ct^/ electron cre- 
ation operators in the SWNT and in lead I, respec- 
tively, and T/(r) describes the transparency of the tun- 
neling contact I. Finally, i/gate — —jJ-gJ^c accounts for a 
gate voltage capacitively coupled to the SWNT, with Mc 
counting the total electron number in the SWNT. 
SWNT Hamiltonian. In the following the focus is on arm- 
chair SWNTs at low energies. Then, if periodic bound- 
ary conditions are applied, only the gapless energy sub- 
bands nearby the Fermi points F = ±Kq = ^KqCx with 
ix along the nanotube axis, are relevant [H, |^. To each 
Fermi point two different branches r = R/L are asso- 
ciated to the Bloch waves (P/j/L fk(^) ~ s*'''^'Pi?/L,F(^)j 
where k measures the distance from the Fermi points 
i^o 13 (Fig. la left). In this Letter, however, we are 
interested in finite size effects. Generalizing to the 
case of SWNTs we introduce standing waves which fulfill 
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open boundary conditions (Fig. la right): 



R/L,K 



V2 



(3) 

with quantization condition k = 'K^m^ + l\)/L^ an 
integer, and L the SWNT length. The offset parameter 
A occurs if Kq ^ Trn/i, and is responsible for the energy 
mismatch between the R and L branches. Including the 
spin degree of freedom, the electron operator reads 



f=R.,L "■'^ 



OBC I 



(4) 



with Cf<jK the operator which annihilates [(ySfP^^*^) \(t). 
The interacting SWNT Hamiltonian then reads 

Hq = hvF^SgTL{f)^K4.^„CfKa + (5) 

with vf the Fermi velocity. We introduce ID operators 



2L ^ 



in terms of which the electron operator in becomes 

*^(^) = E ^SniP)y^fsgn(F}r,F{r)lljfFaix), (6) 
F,f=± 

where we used the convention that R/L = ±1, R/L = 
±1. Upon inserting © into (0, integration over the co- 
ordinates perpendicular to the tube axis yields the inter- 
acting Hamiltonian expressed in terms of ID operators 
and an effective ID interaction K//(a;,a;')- Using stan- 
dard bosonization techniques 0, y| Hq can now be diag- 
onalized when keeping only forward scattering processes 
associated to Veff{x,x'). It reads 





FIG. 1: Energy spectrum of a SWNT with open boundary 
conditions (right) described in terms of left (L) and right 
(7?) branches. It is constructed from suitable combinations 
of travelling waves whose spectrum is shown on the left side. 
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Asgn(f)7Vf, 



E E E "^jSqals^ajSq 

q>0 j—c,s (5— ± 



(7) 



where the first line is the fermionic contribution and rep- 
resents the energy cost, due to Pauli's principle and the 
Coulomb interaction, of adding new electrons to the sys- 
tem. Specifically, Afpa = ^n'^'a-K'^rcrK, is the operator 
that counts the number of electrons in the (fa)-branch, 
Nc = J2fa-^fa yields the total electron number, and 
Eq — hvFY free-particle level spacing (see Fig. 1). 

The term Ec = Wqq is the SWNT charging energy, where 
Wqq — -p- Jq dx dx'Veff{x,x')cos{qx)cos{qx'). The 
second line of Q describes bosonic excitations in terms 
of the bosonic operators ajsq . Four channels are associ- 
ated to total {jS = c+, s+) and relative {j6 = c— ,s— ) 
(with respect to the occupation of the R and L branch) 
charge and spin excitations. Generalized spin-charge sep- 
aration occurs, since for three of the channels the energy 
dispersion is the same as for the noninteracting system, 
£jSq = hvFQ = hvFj^nq =: SQUq, [uq & positive integer), 
but the {c+) channel is affected by the interaction with 
Ec+q — £o"-g(l + ^Wqq/ 80)^^"^ . The cigenstatcs are 



n (^iSq\y 

q>0,jS 



1/2 



iV,0 



N,mj, (8) 



where |A^, 0) has no bosonic excitations and N defines 
the number of electrons in each of the branches (fa). 

Dynamics. Our starting point to describe transport in 
SWNTs is the exact equation of motion 

= Tri^ads [Hi.{t),W\t)] , (9) 

for the reduced density matrix (RDM) — Tri^adsW^ 
of the SWNT. Here {t) is the density matrix of the 
whole system consisting of the leads and the quantum 
dot, and Tr leads indicates the trace over the lead de- 
grees of freedom. The apex / denotes the interaction 
representation with Ht from (j^J as the perturbation. 
We make the following approximations: i) We assume 
weak coupling to the leads, and treat Ht up to sec- 
ond order, i.e., we consider the leads as reservoirs which 
stay in thermal equilibrium and make the factorization 
ansatz {t) = p'{t)psPd =■ p'{t)pieads where p^j^ = 
Z~'l'^e~^^^^i''^^^i''^^i''\ with Zsjd the partition function 
and /3 the inverse temperature, ii) Being interested in 
long time properties, we can make the so called Markov 
approximation, where the time evolution of p^{t) is only 
local in time, iii) Since we know the eigenstates |A^, m) 
of it is convenient to calculate the time evolution 
of p^ in this basis. We assume that matrix elements be- 
tween states representing a different number of electrons 
(charge states) in the SWNT and with different energies 
vanish. Coherences between degenerate states with the 
same energy E are retained! Hence we can divide p^(t) 
into block matrices pll^^ {t\ where i?, N are the energy 
and number of particles in the degenerate eigenstates 



3 



|m). We arrive at equations of the Bloch-Redfield form ergy levels in lead I, and /(e) the Fermi function. Alike, 



kk' 



E E E^fmf.V/..^^MO, (10) 



E' M=N±l kk' 



fa 

/ /"CXj 



where fc, A;' run over all degenerate states with fixed par- 
ticle number. The Redfield tensors are given by {I — s, d) 



R 



kk' - 

I E',MJ 



jjm 



(11) 

and i?f^^fi/ = E,.=±rJJ,^^,^", where the quantities 
p(a)_Biv s„ g^^g transition rates from a state with to 
a state with M particles. Known the stationary density 
matrix pl^, the current (through lead /) follows from 



J — ZXVt, ^ \^ l,njjk J ykn,st■ 

N,E,E' nkj 

(12) 

iv) We exploit the localized character of the transparen- 
cies Ti{r) in Eq. jSJ, and make use of the slowly vary- 
ing nature of the operator ipfFa{x) in Eq. ((HJ. This 
enables us to evaluate the ID operator at the SWNT 
contacts and pull it out from the space integrals which 
enter the definition of the transition rates. It holds 
{r\',l^,„F{x = 0)|s) := (^,.)f/''"+^; (r|Vf<.F(x = L)\s) = 
g-i7rsg„(_F){A^,„sg„(f)+A} (^-^)^j^'^«+i foj. tijg matrix ele- 
ments between the states |r), |s) with energy E' and 
particle number N, -I- 1, respectively. We thus can 
introduce 

^irr'is)^ I d\ I d\'Ti{r)Ti{r')Y^4>i^{TWi^{T') 

X E Sgn(Fi^')<^sgn(i=^)f.F(?^<Psgn(F')f',F'(''')'7i(A): 
FF' 



to describe the influence of the geometry of a tun- 
neling contact at the tube end. The term ?//(A) = 
gi7rsgn(F-F')A(i-5,,rf) accouuts for the mismatch A. As- 
suming a 3D electron gas in the leads, e.g. of gold, we find 
that for a realistic range of energies is <I>;ff' (e) = Sffi^i^ 
i.e. the leads are " unpolarized" . We thus obtain 



^(±)Fjv E'fj + i _ 1 



" Lrss'r' 



^lall depf{e)<^i{e) (Ac 



\En Sjv + i 



with pf{e) — pi{e)f{e), where pi{e) is the density of i 



with pp(£)=p,(£)(l- /(e)). 

When are coherences needed? Eqs. pU|l with p2l) show 
that coherences (in the energy basis) enter the evalua- 
tion of the current. In the low bias and temperature 
regime kBT,eV := e(T4 — Vd) ^ eo, however, where 
only ground states contribute to the current, because of 
(A^,0 \4a^fa\N',0) = [l/2L)5j^,j^, only diagonal ele- 
ments of the RDM contribute. Hence, due to the "unpo- 
larized" character of the leads, the commonly used mas- 
ter equation (CME) with population's dynamics only is 
valid. At larger biases coherences should be included Il4|. 

In the following we focus on the case A w 0, relevant 
to explain the experimental results for sample C in joj. 

Low bias regime (CME is valid). At low bias the cur- 
rent can be obtained by looking to transitions between 
ground states with A^ and A^ -I- 1 particles and energies 
E'^,E^j^^. Then, the matrix element {ipfa)^}^^^ is non 
zero only if A^' = A^ — Cra , with Cra the unit vector, and 



EE(^^-)wiv'(^'=-)jv'w 



2L 



Cn. 



N' 



(15) 



N' 



Here Cm.N' is the number of ground states with A^' par- 
ticles whose configurations A^' differ from the fermionic 
configuration of a given ground state with A^ electrons 
only by a unit vector. With A^ = 4m, 4m -|- 1,4m -|- 
2,4m -I- 3 one finds Cn,n+i — 4,3,2,1 and Cn+i,n = 
1,2,3,4, respectively. We also notice that all ground 
states with A^ particles are populated with equal proba- 
bility, such that we can introduce the occupation prob- 
ability PN{t) = d^P]^ 7v('')' ■where d^ is the degeneracy 
of the ground states with particle number A. It holds 
rfw/rfiv+i = Cn+i,n / Cn,n+i- The corresponding CME 
for PnH) can now be easily solved and the current eval- 
uated in analytic form. We find 



ln,n+i\ 



eA/ CN,N+lCN+l,Nlsld 



J2i=s d li[CN,N+if{£i) + Cn+i.n{1 - /(£;))] ' 

(16) 

where A/ = |/ (e^) - / (ed)|, e; = eVi - AE, and 
AE = E% - E%_^^. Moreover, 7, = {n/ Lh)^ipi{0). 
This expression can be further simplified in the regime 
\eV\ ^ kT <C eo where the linear conductance Gn,n+i 
is obtained by linearizing A/ in V , and by evaluating the 
remaining quantities in (|16|l at zero bias. The conduc- 
tance trace exhibits four-electron periodicity (Fig. 2a), 
with two equal in height central peaks for the transitions 
]\[ = 4rn+l^N + l, N = 'im + 2^N+1, and two 
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smaller peaks for N = Am ^ iV+1, N = 4m + 3 N + 1 
also equal in height. The relative height between central 
and outer peaks is Gf-l,^ im+2/GT^Am+i = 27/(10 + 
4-^6) « 1.36, independent of the ratio js/ld- 

In the bias regime Eq > \eVi ± AE\ > fcT is A/ = 1. 
If e.g. eVs-AE < and eVd-AE > 0, such that tunnel- 
ing is preferable rom source to drain, we find /jv,jv+i = 
eCN, n+iCn+i, Nisid/ {1sCn,n+i + IdCN+iM)- In this 
regime, the nonlinear conductance will still exhibit four- 
electron periodicity, Fig. 2b. For 7^ = jd one still expects 
two equal central peaks and two smaller outer peaks with 
ratio 3/2. li is ^ Id this latter symmetry is lost. If we 
invert the sign of the bias voltage, the current is obtained 
by exchanging 7^ with 7^. 
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FIG. 2; a) Conductance vs. electrochemical potential in the 
linear regime eV <^ ksT <^ eq. Despite asymmetric contacts, 
the two central peaks and the two outer peaks have equal 
height, b) Current in the regime eo S> \eVi\ 3> kT. Asym- 
metry effects become visible. Four-electron periodicity is still 
observed. Parameters are Ec — 9.5 meV; ksT = 0.10 meV, 
£0 = 2.9 meV, 7^ = 5jd = 4.9 ■ lO^^s-V 

High bias regime. In the bias regime eV > Eq states 
with bosonic as well as fermionic excitations contribute 
to transport. An analytical treatment is not possible, 
except to define the position of the various excitation 
lines. The resonance condition for tunneling in/out of 
lead I is as usually given by eVi + AE^ = 0, where 
AE^ = ± (£;Ar±i - En) ■ Besides the resonance condi- 
tion, also the overlap integral between initial and final 
state determines the rates, and hence the "active" reso- 
nance lines contributing to the current. Fig 3a shows the 
current in a bias voltage-electrochemical potential plane 
for the symmetric case 7s = Id- By choosing the addi- 
tion energies provided in |^ A/ii — Ec — 9.5 meV and 
A/i2 = 13.4 meV, we can reproduce all the excitation 
lines from sample C in 0. Moreover, we find a level 
spacing eo = 2.9 meV, which well agrees with the esti- 
mated length for sample C of 750 nm. We compare with 
the mean field parameters: to fit the data, an unusually 
large exchange interaction J ~ 2.9 meV as well as a band 
shift EqA « J had to be assumed in |^ (in our theory 
is EqA « J « 0). This yields a level spacing three times 



larger than the one obtained from our treatment and not 
consistent with the measured SWNT length. 

Finally, the effect of the coherences induced by the 
bosonic excitations is shown in Fig. 3b, where a differ- 
ence plot for the current with and without coherences 
is shown. Though the coherences do not qualitatively 
change the current, they do have a quantitative influence 
in a region of intermediate bias V. A further indication 
for non- Fermi liquid behavior could lie in negative differ- 
ential (NDC) features originating from spin-charge sepa- 
ration, as predicted for a spinful Luttinger liquid quan- 
tum dot [l5j. Asymmetric contacts are a necessary re- 
quirement. We checked these predictions as possible ex- 
planation of the NDC seen in joj. We confirm that (also 
for non-relaxed bosons) NDC occurs. However, very large 
asymmetries must be assumed. Moreover, in contrast to 
the experiments, all the NDC lines have the same slope. 

To conclude, we discussed linear and nonlinear trans- 
port in SWNT quantum dots using a bosonization ap- 
proach. Our results are in quantitative agreement with 
experimental findings in j^] . Further work to explain the 
nature of the NDC seen in loll is needed. 
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FIG. 3: a) Current in a bias voltage - electrochemical poten- 
tial plane for the symmetric contacts case, b) Difference plot 
of the current with and without coherences. Here fcsT — 0.01 
meV and eo/sc+q ~ 0.21. Other parameters are as in Fig. 2a. 
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